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Abstract 

Coupled-mode systems are used in physical literature to simplify the nonlinear Maxwell and 
Gross-Pitaevskii equations with a small periodic potential and to approximate localized solutions 
called gap solitons by analytical expressions involving hyperbolic functions. We justify the use of the 
one-dimensional stationary coupled-mode system for a relevant elliptic problem by employing the 
method of Lyapunov-Schmidt reductions in Fourier space. In particular, existence of periodic/anti- 
periodic and decaying solutions is proved and the error terms are controlled in suitable norms. The 
use of multi-dimensional stationary coupled-mode systems is justified for analysis of bifurcations of 
periodic/anti-periodic solutions in a small multi-dimensional periodic potential. 

1 Introduction 

Gap solitons are localized stationary solutions of nonlinear elliptic problems existing in the spectral gaps 
of the Schrodinger operator associated with a periodic potential. In particular, gap solitons have been 
considered in two problems of modern mathematical physics, the complex- valued Maxwell equation 

and the Gross-Pitaevskii equation 

iE t = -V 2 E + V(x)E + a\E\ 2 E, (1.2) 

where E(x,t) : R N x R i-> C and V 2 = d 2 + ... + d 2 N . For applications of the complex- valued 
Maxwell equation (jl.ip . E is a complex amplitude of the electric field vector, c is the speed of light, 
and n(x,\E\ 2 ) is the refractive index. The scalar equation (jl.ip is valid in the space of one and two 
dimensions (N =1,2) for so-called TE modes but not in the space of three dimensions (N = 3), where 
the system of Maxwell equations for a vector- valued function E must be used [I7j. For applications 
of the Gross-Pitaevskii equation (|1.2p . E is the mean-field amplitude, a is the scattering length, and 
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V(x) is the trapping potential. The scalar equation (|1.2|) is the mean-field model valid in the space of 
three dimensions (N = 3) and it can be used in the space of one and two dimensions (N = 1, 2) under 
additional assumptions [IB] . 

Stationary solutions of the Maxwell and GP equations are found from the elliptic problem 

X7 2 U + u 2 U + eW(x)U = a\U\ 2 U, (1.3) 

where U(x) : R N i— > C and (u! 2 ,e,a) are parameters. The elliptic problem (|1.3p is related to the GP 
equation (jl.2p by an exact reduction E = U{x)e~ tU) t , where the potential V{x) is represented by 
V = — eW(x). The same problem is related to the Maxwell equation (jl.ip with the refractive index 
n(x, \E\ 2 ) = ng(l + fj,W(x) + z^|i£| 2 ) by the exact reduction E = U(x)e~ lcu)t / n ° , where parameters are 
related by e = uj 2 /j, and a = —u) 2 v. 

Let us consider the elliptic problem (jl.3p with a real- valued bounded potential W(x), which is 
periodic in each variable Xj, Vj. The associated Schrodinger operator L = — V 2 — eW(x) is defined 
on C^°(M. N ) and is extended to a self-adjoint operator which maps continuously H 2 (R N ) to L 2 (R N ). 
Therefore, the spectrum cr(L) is real. Suppose that the absolutely continuous part of the spectrum of 
L has a gap of finite size on the real spectral axis. For parameters inside the spectral gap, localized 
solutions of the elliptic problem (jl.3p were proved to exist in the relevant variational problem [19] . 
(Earlier works on bifurcations of gap solitons can be found in [2lllO[[T8].) According to Theorem 1.1 of 
[19] . there exists a weak solution U{x) in ff 1 ^), which is (i) real- valued, (ii) continuous on x G R N 
and (iii) decays exponentially as |x| — > oo. 

We investigate more precise information on properties of the gap soliton U (x) by working with the 
asymptotic limit of small e. When e = 0, the purely continuous spectrum of Lq = —V 2 is non-negative 
and no finite gaps of ct(Lq) exist. However, when N = 1 and / e C 1, narrow gaps of ct(L) diverge 
from a sequence of resonant points on a real axis and gap solitons may bifurcate inside these narrow 
gaps. The coupled-mode system has been used in the physical literature since the 1980s to characterize 
this symmetry-breaking bifurcation of the spectrum c(L) and to approximate one-dimensional gap 
solitons of the problem (|1.3p with N = 1 [22]. Therefore, our work deals mainly with the case N = 1. 
We consider the potential function W(x) according to the following assumption. 



Assumption 1 Let W(x) be a smooth 2tt -periodic function with zero mean and symmetry W(— x) = 
W{x) on x G R. The Fourier series representation ofW(x) is 

W (x) = ^2 w 2m e imx , such that ^ (1 + m 2 ) s \w 2m \ 2 < oo, Vs > -, (1.4) 
where wq = and W2m = W-2m = W2mj Vm £ N. 



It will be clear from Proposition [2] that the sequence of resonant points of o~(L) for L = —d 2 — eW(x) 
is located at u = tu n = ^, n G Z, so that a small gap in the spectrum cr(L) bifurcates generally from 
each point uo = u n , n G N and a semi-infinite gap exists near to = ujq = 0. A formal asymptotic solution 



of the elliptic problem (jl.3p near a resonant point uj 



UJr, 



n G N, is given by 



U{x) = yje a{€x)e 2 +b(ex)e 2 + 0(e) 



n 



+ efi + 0(e 2 ), 



(1.5) 



2 



where the vector (a, b) : C 2 satisfies the coupled-mode system with parameter O € M 



ind + Qa + tD 2n 6 = <x(|a| 2 + 2|6| 2 )a, 
-ink/ + Qb + w 2n a = a(2\a\ 2 + \b\ 2 )b 



2 ,|;,|2 w, I 1 " 6 ) 



and the derivatives are taken with respect to y = ex. 

The coupled-mode system (|1 .6j) can be used for analysis of bifurcations of periodic and anti-periodic 
solutions near narrow gaps in the spectrum cr(L). When the solutions (a, b) of system (|1.6p are y- 
independent and the representation (jl.5p is used, the solution 4>(x) of the elliptic system (jl.3p at the 
leading order is periodic in x for even n £ N and anti-periodic for odd n E N. It follows from system 
(jl.6p that a general family of y-independent small solutions (a, b) near the point (0, 0) has the form 
a = ce i8 and b = ±ce tS , where c £ 1, 6 K, and the nonlinear dispersion relation holds in the form 

n ± w 2n ~ 3crc 2 = 0. (1.7) 

The dispersion relation (jl.7p with c = shows that the eigenvalues for periodic/anti-periodic solutions 
in the linear spectrum of the operator L = —d 2 + eW(x) diverge from the point u 2 = ^ and result in 
a narrow gap in the spectrum o~(L) which lies in the interval 

n 2 2 n 2 

-r ~ \ew 2 n\ < w <Y + e l eu; 2n|- (1.8) 

The nonlinear dispersion relation (jl.7|) with c / shows that the nonlinear periodic/anti-periodic 
solutions bifurcate to the right of the boundaries f2 = ±W2 n for o~ = +1 and to the left for a = — 1. We 
justify the persistence of the leading-order results fll.T|> and (jl.8p by using the method of Lyapunov- 
Schmidt reductions in the discrete weighted space ^ 2 (Z') equipped with the norm 

1 + ^-) \U m \ 2 <oo. (1.9) 

Here Z' is a set of either even or odd numbers and U is a set of Fourier coefficients {f7 m }meZ' i n the 
Fourier series 

I7(z) = yfe V C/ m e>^, C/ m = -1= / W C/(x)e-^ m ^x, (1.10) 

where the factor y/e is introduced for convenience. The function U (x) is periodic if the set Z' is even 
and it is anti-periodic if U is odd. By the Sobolev inequality, if U £ £ 2 (Z') with s > 5, then the 
series (ll.lOp converges absolutely and uniformly in the space of bounded continuous functions C®(M) 
according to the bound 

where we have used the inequality |a||6| < |a| 2 + jH 2 - Our main result on bifurcations of periodic/anti- 
periodic solutions is summarized below. 

Theorem 1 Let Assumption^ be satisfied. Fix n € N ; such that W2 n 7^ 0. Let lo 2 = ^ + eVL, where 
Q, £ R. TTie nonlinear elliptic problem U.3\) with N = 1 /10s a non-trivial solution U(x) in the form 
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hi. lb]) with U G IgC^') f or an V s > 5 anc ^ sufficiently small e if and only if there exists a non-trivial 
solution for (a, b) G C 2 of the bifurcation equations 



Qa + w 2n b - o-(\a\ 2 + 2|6| 2 )a = eA e (a, b), 
06 + w_ 2n a - a(2\a\ 2 + \b\ 2 )b = eB e (a, b), 



2 , ,^_, D ,„ m (1-12) 



where A e {a,b) and B e (a,b) are analytic functions of e near e = satisfying the bounds 

V|e| < e , V|a| + |b| < 5 : |A(a, 6)| < C A (\a\ + \b\), \B e (a, b)\ < C B (\a\ + |6|), (1.13) 

/or sufficiently small eo > 0, /ized 5 > 0, and some constants Ca,Cb > which are independent of e 
and depend on 5. Moreover, A € (a,b) = B e (b,a), V(a,o) G C 2 and 

V|e| < e : U(x) - y/i (ae" + be— J ^ < Ce 3 / 2 . (1.14) 
/or some e-independent constant C > 0. 

Corollary 1 T/ie coupled system admits a symmetry reduction a = b, where the value o/oeC 

satisfies the scalar equation 

Oa + n?2nO — 3cr|a| 2 a = e^4 e (a, a). (1-15) 
Under the reduction, the solution U{x) is real-valued. 

The results of Theorem Q] and Corollary [T] justify the use of the y-independent coupled-mode system 
(|1.6|) for bifurcations of periodic/anti-periodic solutions of the nonlinear elliptic problem (|1.3j) with 
AT = 1. In particular, the only non-zero solutions of the scalar equation (|1.15p occur for either a£l 
or a G iM., when the scalar equation (|1.15|) is reduced to the extended nonlinear dispersion relation 

n±w 2n -3o-\a\ 2 = eA±(\a\), (1.16) 

where A-)-(|a|) = ^A £ (a,a) is a bounded, real-valued error term for sufficiently small e and finite value 
of \a\ G R. Note that the values of A± are real- valued due to the gauge invariance of the coupled-mode 
system (|l,12p inherited from the gauge invariance of the elliptic problem (jl.3p . 

The newly formed gap (jl.8p of the continuous spectrum of L = —d 2 + eW{x) corresponds to the 
interval |fi| < |u>2n|- For instance, let a = +1 and W2 n > 0, then the localized solution of the coupled- 
mode system (jl.6p can be written in the exact form [6j [22] 

a(y) = b(y) = ^ - /< " f (1.17) 

A/3 yW2n - OcOsh(KJ/) + iy/W2n + O Smh(K£/) 

where k = ^v^in ~~ The exact localized solution can easily be found for a = — 1 and W2 n < 0. 
The trivial parameters of translations of solutions in y and arg(a) are set to zero in the explicit solution 
(|1.17p . such that the functions a(y) and b(y) satisfy the constraints a(y) = a(—y) = b(y). 

Definition 1 The gap soliton of the coupled-mode system il.6\) is said to be a reversible homoclinic 
orbit if it decays to zero at infinity and satisfies the constraints a(y) = a(—y) = b(y). 
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We justify the persistence of the gap soliton (|1.17|) of the coupled- mode system (|1.6p in the nonlinear 
elliptic problem (jl.3p by working with the Fourier transform of U(x) 



U(x) = [ U(k)e lkx dk, U(k) = I U(x)e~ ikx dx, (1.18) 

V27T Jr v27re 



where again the factor y/e is introduced for convenience. We develop the method of Lyapunov— Schmidt 
reductions in the continuous weighted space L^(M) equipped with the norm 

UeLl(R): \\U\\ L i m = I (1 + k 2 ) q/2 \U(k)\dk < oo. (1.19) 

By the Riemann-Lebesgue Lemma, if U G £g(R) for some q > 0, then the n-th derivative of the function 
U(x) is bounded and continuous for < n < [q] and it decays to zero at infinity, i.e., U G C^(M) and 
lim U^(x) = 0. Indeed, for any G(x) = U^(x) and G G L X (R), it follows that 

| ac | • — >oo 

\\G(x)\\ c o m <C\\G(k)\\ Lim , (1.20) 

for some C > 0. In addition, the Schwartz space is dense in L 1 (R) such that there is a sequence 
{Gj(k)}j e ^ in the Schwartz space which converges to G(k) in L 1 -norm, and therefore, there exists a 
sequence {Gj(x)}j £ N which converges to G(x) in C^(lR)-norm, such that lim G(x) = 0. 

\x\— >oo 

Related to the coupled-mode system for (a, b) in variable y, we shall also use the Fourier transform 
for (d, b) in variable p, where y = ex and p =\- We note that the norm L 1 (R) is invariant as follows: 

A(k) = -a ( -] : ||^|Ui(«) = ||a|Ui(RV (1.21) 



This invariance explains the choice of the space Lg(M) in our analysis (see also [H]). Our main result 
on the existence of gap soliton solutions is summarized below. 



Theorem 2 Let Assumption [1] be satisfied. Fix n G N, such that ui2 n 7^ 0. Let u = %- + eft, such 
that < |i02n|- Let a(y) = b(y) be a reversible homoclinic orbit of the coupled-mode system U.6\) in 
Definition^ Then, the nonlinear elliptic problem ll.3\) with N = 1 has a non-trivial solution U(x) in 
the form il.18]) with U G L x (R) for any q > and sufficiently small e such that 

V|e| < e : \\U{x)-yfe 

for some sufficient small cq and e-independent constant C > 0. Moreover, the solution U{x) is real- 
valued, continuous on x G R, and lim U(x) = 0. 

|z|^oo 



mx 

a{ex)e 2 + a(ex)e 



2 



(1.22) 



The results of Theorem [2] give more precise information about gap solitons of the elliptic problem 
(|1.3p with N = 1 compared to the general result in Theorem 1.1 of [19J, since the leading-order 
approximation of U(x) is given by the exponentially decaying solutions (|1.17p of the coupled-mode 
system (ll.6p . On the other hand, we do not prove in Theorem [2] that U(x) decays exponentially as 
\x\ — > 00. 
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Rigorous justification of the approximation (|1.5|) and the time-dependent extensions of the coupled- 
mode system (|1.6p were developed in [12! for the system of cubic Maxwell equations and in [21j for the 
Klein-Fock equation with quadratic nonlinearity. A bound on the error terms was found in the Sobolev 
space H X (R) in [12] and in the space of bounded continuous functions CP(R) in [21]. The bound is 
valid on a finite interval of the time evolution, which depends on e. The error is not controlled on the 
entire time interval t E K and, in particular, the formalism cannot be used for a proof of persistence of 
the leading-order approximation (jl.5p and (|1.17p for the stationary solutions of the nonlinear elliptic 
problem f)l .3|) . The results of our Theorem[2]are more precise than Theorem 1 of [12] and Theorem 2.1 of 
[21] in this sense, since the error bound of the leading-order approximation is controlled independently 
of t E K. In a similar context, the justification of the nonlinear Schrodinger equation for the nonlinear 
Klein-Gordon equation with spatially periodic coefficients is reported in [4]. 

The method of Lyapunov-Schmidt reductions for periodic solutions was used in [8j[IT]. The work [8] 
deals with a two-dimensional lattice equation for the nonlinear wave equation when eigenvalues of the 
relevant linearized operator accumulate near the origin. In this case, the Nash-Moser Theorem must 
be used for the infinite-dimensional part of the Lyapunov-Schmidt decomposition. In our case, the 
linearized operator for the one-dimensional lattice equation has eigenvalues bounded away of the origin 
and the Implicit Function Theorem can be applied without additional complications. This application 
of the technique is similar to the one in [11], which deals with the periodic wave solutions in the system 
of coupled discrete lattice equations. Other applications of the method for periodic wave solutions in 
equations of fluid dynamics can be found in [51 [7] . 

Persistence of modulated pulse solutions was considered in [131 E] i n the context of the nonlinear 
Klein-Gordon equations. (Earlier results on the same topics can be found in [3j [9].) Methods of spatial 
dynamics were applied to a relevant PDE problem for modulated pulse solutions, the linearization of 
which possessed infinitely many eigenvalues on the imaginary axis. The local center-stable manifold 
was constructed for the nonlinear Klein-Gordon equations after normal-form transformations and the 
pulse solutions were proved to be localized along a finite spatial scale, while small oscillatory tails 
occur generally beyond this spatial scale. In contrast to these works, we will not reformulate the ODE 
problem as an extended PDE problem and avoid the construction of the local center-stable manifold. 
This simplification is only possible if the variables of the time-dependent problems (11.11) and (II. 2p can 
be separated and modulated pulse solutions are described by the reduction to the elliptic problem (jl.3p . 
We note that the basic equations of electrodynamics, such as the real- valued Maxwell equation (of the 
Klein-Gordon type) , would not support the separation of variables and the modulated pulse solutions 
do not generally exist in the real- valued Maxwell equation [3]. 

The article is structured as follows. The proof of Theorem [1] is given in Section 2, where the 
technique of Lyapunov-Schmidt reductions in ^(Z') is developed for bifurcations of periodic/anti- 
periodic solutions. The proof of Theorem [2] is given in Section 3, where the technique of Lyapunov- 
Schmidt reductions in is extended for persistence of decaying solutions. Section 4 discusses 
applications of similar methods for the justification of multi-dimensional multi-component coupled- 
mode systems with N > 2. 

2 Lyapunov-Schmidt reductions for periodic/anti-periodic solutions 

Let the potential W{x) and the solution U(x) to the elliptic problem (II. 3p with N = 1 be expanded 
in the Fourier series (|1.4j) and (|1.10p respectively. By using the Fourier series, we convert the elliptic 
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problem (|1.3|) with N = 1 in the space of bounded, continuous, periodic/anti-periodic solutions U G 
C°(1R) to a system of nonlinear difference equations in the discrete Sobolev weighted space U G £ 2 (Z') 
for some s > ^. The nonlinear difference equations are written in the explicit form 



fw 2 —jU m + e ^2 w m - mi U mi = ea ^ ^ f7 mi J7_ m2 C/" OT _ mi _ m2 , Vm G Z', (2.1) 

which can be casted in the equivalent matrix-vector form 

(£ + eW)U = ecrN(U,U,U). (2.2) 

Here U is an element of the infinite-dimensional vector space Zf(Z') with the norm (|1.9p . elements of 
matrix operators £ and W are given by 

(TTZ 2 \ 
- — J 5 m ,fc, W mifc = f m -t, V(m, fc) G Z' x Z', 

and N(U, U, U) = U * 7£U * U consists of the convolution operator with the elements (U * V) m = 

Yl UkV m ^k and the inversion operator with the elements (TZXJ) m = U- m . 

fcez 

We shall verify that the nonlinear vector field associated with the difference equations (|2.ip is closed 
in space U G ^ 2 (Z') with s > ^ (Lemma [T]). Working in this space, we shall apply the Implicit Function 
Theorem in two cases u ^ M\{^} ne z and uj = uj n = ^ for some n G N. The first case is non-resonant 
and the Implicit Function Theorem guarantees existence of the unique zero solution U = of system 
(12. 2|) near U = and e = (Lemma [2]) . The second case corresponds to a bifurcation of non-zero 
periodic or anti-periodic solutions of system (|2.2p and it is analyzed by using the Lyapunov-Schmidt 
decomposition. To prove Theorem [TJ we will prove that there exists a unique smooth map from the 
components (U n , U- n ) to the other components U m , Vm G Z'\{n, — n}. Projections to the components 
(U n ,U- n ) yield the coupled-mode equations (II. 121) . while the bounds on the remainder terms (|1.13p 
follow from the bounds on the vector U in space l 2 (TJ). Representation (|1.14p and symmetry reductions 
of Corollary [T] follow from the technique of Lyapunov-Schmidt reductions. 



Lemma 1 Let W G £ 2 (Z) for all s > \. The vector fields W* and N map elements of / 2 (Z') with 
s > i to elements ofl 2 (Z'). 

Proof. The space Z 2 (Z) with s > ^ forms a Banach algebra. Therefore, there exists a constant 
< C(s) < oo such that 

VU,VG/ 2 (Z): ||U*V||,; (z) <C7( a )||U|| e(Z j||V||,; (z) , Vs > \. (2.3) 

This property was proven in [8] and it is similar to the one in the continuous H S (M) spaces. The maps 
W * U and U * TZ\] * U act on U G Z 2 (Z'), where Z' is a set of either even or odd numbers. Both 
convolution operators transform a vector on Z' to a vector on Z'. Therefore, ^ 2 (Z') forms a linear 
subspace in the vector space Z 2 (Z) with the same algebra property (|2.3p . As a result, both W*U and 
N(U, U, U) map elements of Z 2 (Z') with s > \ to elements of Z 2 (Z'). □ 

Lemma 2 Let W G / 2 (Z) with s > h and uj G K\{f} n ez f° r n £ Z. The nonlinear lattice system 
\2.2\i has a unique trivial solution U = in a local neighborhood of U = and e = 0. 
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Proof. If u G R\{^} n6 z, the operator C in system (12. 2|) is invertible and 

WC^Wq^q < 1 -2- = /> <OO, Vs>0. 

mm \uo 2 - ^jr\ 

It follows from the estimate (II. lip that Yl \ w m-k\ < 00 f° r an y "i€Z' and s > ^. Therefore, the 

matrix operator W is a relatively compact perturbation to C. By the perturbation theory [15], there 
exists an e-independent constant < p < oo such that 

|| (£ + eW)" 1 || M < p, Vs > 0, 

for sufficiently small e. By Lemma [U the vector field on the right-hand side of system (|2.2j) is closed 
in Zg(Z') for s > ^. Moreover, it is analytic with respect to U £ ^s(Z') an d e G R. The zero solution 
U = satisfies the nonlinear lattice system (|2.2p for any e G M. The Frechet derivative of system (|2.2p 
at U = (which is just the operator £ + eW) has a continuous bounded inverse for sufficiently small 
e. By the Implicit Function Theorem, the zero solution U = is unique in a local neighborhood of 
U = and e = 0. □ 

Proof of Theorem^ If oj = ^ for some n G N, the operator C is singular with a two-dimensional kernel 

Ker(£) = Span (e n , e_ n ) C Z 2 (Z'), 

where e n is a unit vector in l 2 (U). The straightforward decomposition of / 2 (Z') = Ker(£) © Ker (£)-*- 
is nothing but the representation 

U = ae n + 6e„ n + g, (2.4) 

where 

g G KerOC)^ = {g G Z 2 (Z') : <7n = <?-n = 0}. (2.5) 

Let V be the projection operator from / 2 (Z') to Ker (£)-*- at w 2 = j. It is obvious that V CP is a 
non-singular operator at uj 2 = By using the same argument as in Lemma [21 we obtain that there 
exists an e-independent constant < p < oo, such that 

\\(V(£ + eW)rr 1 \\ l 2^ l 2<p, Vs>0 (2.6) 

for sufficiently small e. The inhomogeneous problem for g is written in the explicit form 

h eSl j g m + e w m ^ k g k 

- ecr U mi U- m2 U m - mi - m2 = - e (a-w m _ n + 6-u; m+n ) , (2.7) 

mieZ' m 2 ez' 

where Z7 m = g m + a<5 m) „ + b5 m - n . By Lemma [U the vector field of system (12. 7|) is closed in / 2 (Z') for 
s > \ and any (a, 6) G C 2 . Moreover, it is analytic with respect to g G Z 2 (Z') for all (a, b) G C 2 and 
e G M. By the bound (|2.6p and the Implicit Function Theorem there exists a unique trivial solution 
g = of system (|2.7p for (a, b) = (0, 0) and any e G K. It is also obvious that the zero solution exists 
for any (a, b) G C 2 and e = 0. For all (a, 6) / 0, the Frechet derivative of system (12. 7h at g = is 
different from the matrix operator V(C + eW)V by the additional terms 

-ecr [(|a| 2 + |&| 2 )<5 m ,fc + ab5 mtk +2n + a&<W-2n] , V(m, fe) G Z' x Z'. 
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For sufficiently small e and finite (a,b) £ C 2 , these terms change slightly the bound p in (|2.6[) . such 
that the Frechet derivative operator of system (|2.7p at g = has a continuous bounded inverse for 
|e| < eo, where eo is sufficiently small. By the Implicit Function Theorem, there exists a unique map 
G £ : C 2 i— ► Ker (£)-*- C 1 2 {U), which is analytic in e with the properties G e (0, 0) = and Go (a, b) = 0. 
Therefore, the map G t admits the Taylor series expansion in e. The first term of the Taylor series is 

G e (a, b) = e [ag a + bg b + a 2 bg c + ab 2 g d ] + 0(e 2 ), 

where non-zero components of vectors g a ,b,c,d in the constrained space (|2.5p are 

i \ 4lL' m _ n 4:W m j,- n , . ^O~0~m,3n / \ 4(T(5 m) _3 n 

\ga)m — 9 9) \gb)m — 9 ni \gc)m — 9 oi \gd)m — 9 9" ■ 

m z — n A m z — n A — n z m — n A 

Let \a\ + \b\ < S and 5 is fixed independently of e. Due to the analyticity of G e in e, there exists an 
e-independent constant C > such that 

V|e|<e : ||G e (a, 6)11^ < eC (|o| + |6|) . (2.8) 

The projection equations to the two-dimensional kernel of C is found from system (|2.1[) at m = ±n in 
the explicit form 

—mi—m,2 

miSZ' m 2 eZ' 

(£1 + wo)b + w-2 n a- a ^ S U mi U- m2 U- n - mi - m2 

miSZ' m 2 eZ' 

where f7 m = a<5 mjn + bd m - n + <? m and the map g = G e (a, 6) is constructed above. At e = 0, we obtain 
that g = and 

(|a| 2 + 2|6| 2 )a, 
(2\a\ 2 + \b\ 2 )b. 

This explicit computation recovers the left-hand side of system (|1.12p . The right-hand side is estimated 
from the bound (|2.8[) on the map G e (a, b) in l 2 (Z') with s > \ to yield the bound (fTT3l) . 

To prove the last assertion of Theorem HJ we shall prove that the map G e (a, b) has the symmetry 
(G t ) m (a,b) = (G e )_ m (6, a). Indeed, since W(x) is real-valued, its Fourier coefficients satisfy the 
constraint W-2 m = W2 m , Vm £ Z. The systems of equations ()2.7p and (12. 9p are symmetric with respect 
to the transformation (a,b,g m ) t— * (b,a,g~ m ). By uniqueness of solutions of system (|2.7p in a local 
neighborhood of e = 0, we obtain that (G e ) m (a, b) = (G e )__ m (6, a), Vm G Z'. Then, it follows directly 
from system (pT9|) that A e (o, 6) = -B e (&, a). □ 

Proof of CorollaryUl Due to the property j4 e (a, b) = -B e (6, a), system (|1.12p has the symmetry reduction 
b = a, which results in the scalar equation (|1.15|) . The vector U is given by the decomposition (12. 4p 
with b = a and (G e ) m (o, o) = (G e )_ m (a,o). Therefore, the solution U(x) recovered from the Fourier 
series (ll.lOp is real-valued. □ 



Wn-k9k, 

fcSZ'\{n,-n} 
>'\{n,—n} 



(2.J 



miSZ' m 2 6Z' 

^ ^ c/ mi C7_ m2 C/_ 

miSZ' m 2 eZ' 



-n—mi — m 2 
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3 Lyapunov— Schmidt reductions for gap solitons 



Let the solution U(x) to the elliptic problem (jl,3p with N = 1 be represented by the Fourier transform 
(|1.18p . while the potential W(x) is given by the Fourier series (11. 4p . The elliptic problem (|1.3p with 
N = 1 is converted to the integral advance-delay equation for the Fourier transform U(k): 



or 



k 2 ) U(k) + e ^ w 2m U(k - m) = ea- 



rn <E i 



U{ki)U{k 2 )U{k -h + k 2 )dk 1 dk 2 



Vk G 



(3.1) 



Working in the Fourier space U G Lg(R), where the vector field of the integral advance-delay equation 
(|3.ip is closed (Lemma [3]), we decompose the solution U(k) into three parts 

U(k) = U + (k)xw + (k) + U-(k)xw_(k) + U (k) XK (k), (3.2) 

where X[a,b](k) is a function of compact support (it is 1 on k G [a,b] and on k G M\[a,b]) and the 
intervals m / 



.. ™_ and M' 



arc 



,2/3 



2/3 



', U R'_ 



(3.3) 



Here tu n = \ is a bifurcation value of u and the components U±(k) represent the largest part of the 
solution U(k) near the resonant values k = iw n , which is approximated by the solution of the coupled- 
mode system f 1 1 . 6 1) in coordinates y = ex in physical space and p = - in Fourier space. The intervals 
surrounding the resonant values k = ±ui n have small length 2e 2 / 3 , where both the constant c = 2 
and the scaling factor r = | are fixed for convenience. In fact, we could generalize all proofs for any 
constant c > and any scaling factor | < r < 1. 

In order to prove Theorem [21 we shall apply the method of Lyapunov-Schmidt reductions in space 
Lg(M.) with q > 0. First, we prove the existence of a unique smooth map from (U + (k), U-(k)) to Uo(k) 
(Lemma H]). The solutions for the components (£/+(&), £/_(&)) are then approximated by the suitable 
(exponentially decaying) solutions (a(p),b(p)) of the coupled-mode system (|1.6p rewritten in Fourier 
space (Lemma [S]). The approximation yields the desired bound (jl.22p for U G Lg(M) with q > 0. The 
reduction to the real-valued solutions U(x) becomes obvious from the decomposition of the Lyapunov- 
Schmidt reduction method. Continuity and decay conditions on U (x) follow by the Riemann-Lebesgue 
Lemma. 



Lemma 3 Let W G Z 2 + „(Z) for all s > | and q > 0. T/ie vector field of the integral equation \3.1 
maps elements of £i(R) toiift q > to elements of L\ 



Proof. The convolution sum in the integral equation (13. ip is closed due to the bound 



VU G Lj(R), VW G / 2 +(? (Z) : 



X! w 2mU(k 



m) 



< II^IIli(m)I|w|| Z i (z) < ||[/|| L i ( 



l w H^ +9 (z) 



for any q > and s > i, where the inequality (jl.lip has been used. The convolution integral is closed 
due to the bound 

yu, V G Ll(K) : 



l7(Jfei)V(ife - fei)dfc] 



< CIICII T.ifm \\V\ 
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for any q > and some C > 0, which occur in the inequality 

l + {k 1 + k 2 ) 2 < C(l + fc?) (l + fc£) 

for all (ifci.Jfea) G M 2 . 



□ 



2 

Lemma 4 Lei Assumption [7] Z>e satisfied and u = \ + efi, where n G N and 11 6 R. There exists a 
unique map U e : L](R' + ) x L*(R'_) h-> Lj(M' ) /or aZZ g > 0, swc/t tfiai E/b(fc) = Z7 e (Z7 + , Z7_) and 



V|e| < e : ||Z7o(&)IUi(R^) < e 1/3 C (ll^+IUjCR^ + H^-IUj(»'_)) 



(3.4) 



where cq is sufficiently small and the constant C > is independent of e and depends on 5 in the bound 
\U4-W rxrui \ + HtZl ||li(r' ) < 5 for a fixed e-independent 5 > 0. 



Proof. We project the integral advance-delay equation (13. ip onto the interval k G M.' : 



+ eil, - k 2 I U (k) + e ^ n^mXR^)^ - m) = (&) 



n 

T 



where U(k) is decomposed by the representation (|3.2p . Since 

,2 



U{ki)U{k 2 )U{k -ki + k 2 )dk 1 dk 2l 



mm 



77 
T 



A: 



> 0,,e 2 / 3 , 



for some C n > 0, the linearized integral equation at e = is invertible such that 

.2 n -1 



k- 



< 



1 



C„e 2 /3- 



(3.5) 



The linearized integral equation on Uo(k) for e 7^ is given on k G Mq by 



l^ + eQ-k 2 ) U Q (k) + e ^ w 2m xw (k)U (k - m) 



(*) 



TJ+(k x )U-(k 2 ) + U+fap-fa^ Z7 (A; + k x + k 2 )dk x dk 2 

-evXmdk) / / U+{k 1 )U + (k2)Uo(k-ki + k 2 )dk 1 dk2 
Jr' + Jr' + 

(k) [ [ U-{k 1 )U-{h 2 )Uo(k-k 1 + h2)dk 1 dk2. 



Recall that e 2/3 > e for sufficiently small e. If W G Z 2 +? (Z) and U± G Lg(K±) for all s > \ and g > 0, 
the convolution sums and integrals are closed by Lemma [3l Fix ||Z7+||li(r' + ) + II^-IIli(r' ) < & f° r some 
e-independent 5. Then, the linearized operator is continuously invertible for all |e| < eo- The integral 
equation is analytic in e and admits a unique trivial solution Uo(k) = if U±(k) = on k G M!± or if 
e = 0. By the Implicit Function Theorem, there exists a unique map U t : L^(W + ) x Lg(M.'_) 1— ► L^(Rg) 
for |e| < eo- The map is analytic in e near e = with the properties Uq(U+, Z7_) = and Z7 e (0, 0) = 0. 
Due to the analyticity of the map U t (U + ,U-) in e and the bound (|3.5[) on the inverse operator, the 
solution Uo(k) = U e (U + , Z7_) satisfies the desired bound (|3.4p . □ 
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Lemma 5 Let Assumption [I] be satisfied. Fix n G N, such that u>2 n ^ 0. Let oj = ^ + efi, smc/i 
i/iai \£l\ < \w2n\- Let a{y) = b(y) be a reversible homoclinic orbit of the coupled-mode system U.6\) in 
Definition^ Then, there exists a solution of the integral equation ( lg.il) . such that Uo(k) = U e (U+,U-) 
is given by Lemma^ and 



V|e| < e : 



U+{k) - -a I — - — 



LUR> ) 



e V e 



LJ(R'_) 



< C^/a, (3.6) 



/or sufficiently small eo > anc? e-independent constants C a ,Cf, > 0. 



Proof. Let us use the scaling invariance (jl.2ip and map the intervals for U±(k) to the normalized 
interval M = [-e _1/3 , e _1/3 ] for 



a(p) = et/4. 



k - u n 



b(p) = ell- 



k + u r , 



(3.7) 



The new functions a(p) and b(p) have a compact support on p G ]Ro> while the norms ||a||z,i(R ) and 

||6||z,i(R ) are equivalent to the norms ||^+||l1(r' + ) an< ^ )■ Using the bound (13. 4|) . we project 

the integral equation (13. ip to the system of two integral equations on p G Mq : 



(fi + t^o - np) a(p) + W2 n b{p) ~ °~ 



a(pi)a(p 2 ) + b{pi)b{p 2 ) a(p-pi + p 2 )dp 1 dp; 



-a 



a(pi)b{p 2 )b(P -Pi+ P2)dpidp 2 = ep 2 a(p) + e 1/3 A € (a, b, U e (a, b)), (3.8) 



(f2 + w + np) 6(p) + w-2 n a(p) - cr 



a(pi)a(p 2 ) + b(Pi)b(P2) b(p -p\+ P2)dpidp 2 



-a 



b(Pi)a(P2)a(p -pi+ P2)dpidp 2 = ep 2 b{p) + e 1/3 B e (a, b, U e {d, &)), (3.9) 



where A t and -B e are computed from (a, 6) and the map U t (a,b) of Lemma [H When the right-hand 
side of system f)3.8j) (j3.9j) is truncated and the integration is extended to p G M, system f)3.8j) (|3.9|) 
becomes the coupled-mode system (|1.6p rewritten after the Fourier transform in y. If ^ and 
1^1 < If^nli the coupled-mode system (|1.6p has a reversible homoclinic orbit a(y) = b{y). The Fourier 
transform d(p) decays exponentially as \p\ — > oo, such that the integrals of the system (I3.8p - (l3.9p on 
p G M\Mo are exponentially small in e. Therefore, they can be moved to the right-hand side of the 
system. In addition, the remainder terms A e and B e are analytic in e and controlled by the bound (13. 4p 
on (a(p),b(p)) in the space L^(Ro) for all q > 0. The linear terms p 2 a(p) and p 2 b(p) are also controlled 
by the bound 



e||p 2 a 



M<^ 1/3 H« 



5o)' 



e\\p 



to) 



< e 1 / 3 



iJ(«o)- 



(3.10) 



Therefore, system ()3.8p - (|3.9p is a perturbed coupled-mode system fjl .6|) in Fourier space with the 
truncation error of the order ©(e 1 / 3 ) measured in space L*(Rq). (Note that the system (I3.8p - (|3.9p is 
not closed on -Li(M) as the terms p 2 a(p) and p 2 b{p) represent a singular second-order perturbation of 
the first-order coupled-mode system.) 
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To prove the persistence of decaying solutions of the coupled- mode system (|1.6p . we note that the 
linearized differential operator associated to the coupled-mode system (|1.6p in the physical space is 
given by a self-adjoint system of 4-by-4 component Dirac operators: 



indy + Wo 

-<ra? 
W2n — 2aab 

—2aab 



—aa 



2 




—indy + Wo 

—2aab 
W2n — 2aab 



(3.11) 



where Wq = $7 — 2er(|a| 2 + \b\ 2 ). By Theorem 4.1 and Corollary 4.2 in [6J, the linearized operator (|3.1ip 
is block-diagonalized into two uncoupled 2-by-2 Dirac operators, each has a one-dimensional kernel. 
The two-dimensional kernel of the linearized operator (|3.1ip is related to the translational symmetries 
in y and arg(a) with the eigenvectors [a'(y), b'{y), a'(y), b'(y)] T and [ia(y),ib(y),—ia{y),—ib{y)] T . The 
zero eigenvalue of the linearized operator (|3.1ip is bounded away from the continuous spectrum and 
other eigenvalues on the real axis [Bj. The extended coupled-mode system given by the system (|3.8p - 
(|3.9D after the Fourier transform is only solvable if the right-hand-side lies in the range of the linearized 
operator ([3. lip . 

The nonlinear elliptic problem (jl.3p with real-valued symmetric potential W{x) = W(—x) = W(x) 
has two symmetries: the gauge invariance U(x) — > e ia U{x) for all a £ R and the reversibility U(x) — > 
U(—x). The new system obtained after the Fourier transform (|1.18p and the decomposition (|3.2p 
inherits both symmetries, such that the extended coupled-mode system is formulated in a constrained 
subspace orthogonal to the kernel of the linearized operator (|3.1ip . As a result, the linearized operator 
is continuously invertible in the constrained subspace of space Lg(Mo) x -^q(^o) f° r a h 9 > 0. Truncation 
of the integral terms introduces a small error in the remainder terms but does not change the symmetries 
of the extended coupled-mode system and does not alter the invertibility of the linearized operator. 
By the Implicit Function Theorem, there exists a unique solution of system (I3.8p ^ (|3.9p for a(p) and 
b(p) on p £ Mo j which is close to the reversible homoclinic orbit of the coupled-mode system (jl.6p in 
L^-norm. □ 

Remark 1 If W2 n > and £1 = u>2n, the exact solution (I1.17P describes an algebraically decaying 
reversible homoclinic orbit of the coupled-mode system (|1.6p . Since the continuous spectrum of the 
linearized operator (|3.1ip touches the zero eigenvalue in this case, persistence of algebraically decaying 
reversible homoclinic orbits can not be proved in Lemma 

Remark 2 The symmetry condition on the potential W(—x) = W(x) in Assumption [1] is important 
for the proof of persistence of homoclinic orbits in the nonlinear elliptic problem (jl.3p since it ensures 
that the set of homoclinic orbits of the nonlinear problem (ll.3p includes the symmetric (reversible) 
homoclinic orbits U(—x) = U(x). Lemma[5]can not be proved if the right-hand-side of the extended 
coupled- mode system is not in the range of the linearized operator (]3.1ip . which may occur when a 
homoclinic orbit of the coupled-mode system (jl.6p is positioned arbitrarily with respect to the general 
potential function W(x). Non-persistence of such homoclinic orbits is usually beyond all orders in the 
asymptotic expansion in powers of e and it is typical that the gap solitons persist at two particular 
points on the period of the potential W(x) (see [20J for details). In our paper, we avoid the beyond- 
all-orders problem by imposing a symmetry condition on W{x) which is sufficient for existence of a 
reversible homoclinic orbit which is centered at the point x = 0. 
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Proof of Theorem^ When the solution U (x) is represented by the Fourier transform U(k), both scaling 
transformations of Lemma H] and [5] are incorporated into the solution, and the bounds (|3.4p and ()3.6j) 
are used, we obtain the bound 



V|e| < e : 



U(k) 



1 



k - n/2 



e \ e 



< Ce 1 / 3 , 



(3.12) 



which implies the desired bound (|1.22p in original physical space. It remains to prove that the solution 
U(x) is real- valued. This property follows from the symmetry of the map Uo(k) = U e (U+,U-) con- 
structed in Lemma U] with respect to the interchange of U+(k — u) n ) and U-(k + uo n ) and the complex 
conjugation. As a result, system f)3.8f) — (j3.9f) has the symmetry reduction a(p) = b(p), which is satis- 
fied by the solution of the truncated system. When the partition f|3.2[) is substituted into the Fourier 
transform (11.180 with the symmetry U+(k — uj n ) = U-(k + uj n ), the resulting solution U(x) is proved 
to be real-valued. □ 



4 Lyapunov— Schmidt reductions in multi-dimensional potentials 

Let us consider the elliptic problem (|1.3I) in the space of two dimensions (N = 2). Let the potential 
VF(x) be periodic in both variables with the same normalized period, such that 

W{ Xi + 2tt,x 2 ) = W(x 1 ,x 2 + 27r) =W{ Xl ,x 2 ), V(xi,i 2 )gM 2 . (4.1) 

We shall justify the use of multi-component coupled-mode systems for the analysis of bifurcations of 
two-dimensional periodic/anti-periodic solutions of the elliptic problem (jl.3p . We use the Fourier series 
for the potential W(x) and the solution U(x): 

W(x) = w 2m e m - x , U{x) = V~e £ U m e>- X , (4.2) 

where m ■ x = rriixi + m 2 x 2 and the sets and U 2 are even or odd if the solution U{x) is periodic or 
anti-periodic in the corresponding variable x\ and x 2 . The elliptic problem (jl.3p with N = 2 transforms 
to a system of nonlinear difference equations, which is similar to system (I2.ip : 

— j U m +6 ^ ^ W m — mi U mi — 6CT ^ ^ ^ ^ U mi U— m 2Um—mi—m27 (^"^) 



for all m 6 Z'j x Z^. The nonlinear lattice system (|4.3p is closed in the space l 2 s {U x x Z' 2 ) with s > 1 
thanks to the Banach algebra property: 

VU, V G l 2 s (Z 2 ) : ||U*V||,2 (za) < C( a )[|U[|,j (za) ||V||,; (z2)) Vs > 1, (4.4) 

for some C(s) > 0. Under the same constraint s > 1, the double Fourier series (|4.2p converges absolutely 
and uniformly in C®(M. 2 ). 



The system (|4.3p takes the same abstract form (|2.2p . where U is an element of the vector space 
l 2 {^'i x 7j' 2 ) with s > 1. An extension of Lemma [2] tells us that no non-trivial solution U of the system 
3|) exists in a local neighborhood of U = and e = unless u = tu n = where n = (ni,n 2 ) £ 1? 
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and 1 77. | = -\Jn\ + n 2 ,. Bifurcations of non-trivial solutions occur only in the resonant case u> = u> n and 
the number of bifurcation equations (leading to the coupled-mode system) is defined by the dimension 
of the resonant set S n in 

S n = {m€ Zi x Z 2 : \m\ 2 = \n\ 2 } . (4.5) 
Here again the set Zj is even/odd if n\ is even/odd and so is the set Z 2 with respect to ni- 

Lemma 6 The set S n admits the following properties: 

(i) < Dim(S n ) < 00. 

(ii) If n = 0, the zero solution m = is unique. 

(Hi) If n = (m,0), then Dim(S' n ) > 2 if n\ is odd and T)im(S n ) > 4 if n\ is even. 

(iv) If n = (m, n2) G N 2 , then Dim(S' n ) > 4 if n\ — ri2 is odd and Dim(S' n ) > 8 if n\ — n2 is even and 
non-zero. 

Proof, (i) follows from the bound \m\ 2 < 00 on the space of integers and from the existence of the 
solution m = n. (ii) is obvious from \m\ 2 = 0. (iii) follows from the existence of particular solutions 
(±7ii,0) and (0,±ni) of \m\ 2 = n 2 (if n\ is odd, the solutions (0,±ni) do not belong to the space Z 2 
of even numbers), (iv) follows from the existence of particular solutions (±ni, ±712) and (±ri2, ±tii) of 
\m\ 2 = n 2 + n?, (if n\ — 712 is odd, the solutions (±ri2, ±"-i) do not belong to the space U 1 x Z' 2 of the 
opposite parities and if n\ = n 2 , the solutions (±n 2 ,±ni) are not different from (±^,±722)). □ 

Proposition 1 Let W 6 / 2 (Z 2 ) for all s > 1 and uj 2 = ^ + eVL for some n G N 2 and Qel. Let 
the set S n be defined by J^.5[ ) with d$ = Dim(5 n ). The nonlinear lattice system {J^.S) has a non-trivial 



solution U G l 2 (^i X Z 2 ) for all s > 1 and sufficiently small e if and only if there exists a non-trivial 
solution for a G C ds of the bifurcation equations 

^ a i m + Yl w m- mi a jmi -a ^2 a im i aj_ m2 a Jm _ mi „ m2 = eA jm>e (a), Vm G S n , (4.6) 

m 1 eS„ m 1 £S n -m2&S' n 

where j m is an index of m in the set S n , the set S' n is a subset of S n , such that m — mi — ni2 G S n , 
A e (a) G C ds depends analytically on e near e = 0. Moreover, there exists constant C > which is 
independent of eo > and depends on 5 > such that 

V|e| < e , V||a|| /1(C d s) < 5 : || A e (a)|| J1(C d s) < C||a|| /1(C d s) , (4.7) 

where eo is sufficiently small and 5 is fixed independently of eo. 

Proof. The proof repeats the proof of Theorem [1] due to the fact that \m\ 2 — \ n\ 2 are bounded away from 
zero for m G Z' x x 7j' 2 \S n . The Lyapunov-Schmidt reductions are performed after the decomposition 
U = J2mes n a i m e m + g, where g G Ker(£) ± . □ 

Example 1 The resonant value ^(0,0) = corresponds to the single-mode bifurcation, like in the one- 
dimensional problem (N = 1). The next resonant value W(i,o) = w (o,i) = \ corresponds to the two-mode 
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bifurcation, which has the same coupled-mode equations as in the one-dimensional problem (N = 1) 
due to the separation of the periodic Fourier series in one variable and the anti-periodic Fourier series in 
the other variable. Finally, the next resonant value = gives the first example of the non-trivial 
four-component coupled-mode equations in the space of two dimensions (N = 2). The coupled-mode 
equations (|4.6p can be rewritten explicitly for the components (01, a 2 , a 3 , 04) which corresponds to the 
Fourier modes for the resonant set Snu = {(1, 1); (—1, —1); (1, —1); (—1, 1)} at the selected order: 

(12 + u>o,o)ai + ^2,2^2 + ^0,203 + ^2,004 = & ((|ai| 2 + 2|a 2 | 2 + 2|a 3 | 2 + 2|a 4 | 2 )ai + 2a 2 a 3 a4) , 

(f2 + wofi)a 2 + W-2-2CL1 + ^-2,o«3 + ^0,-204 = cr ((2|ai| 2 + |a 2 | 2 + 2|a 3 | 2 + 2|a4| 2 )a 2 + 2aia 3 a 4 ) , 

(O + u>o,o)a3 + w^2,-2«4 + 1^0,-201 + ^2,002 = cr ((2|ai| 2 + 2|a 2 | 2 + |a 3 | 2 + 2|a 4 | 2 )a 3 + 2a 4 aia 2 ) , 

(O + u> ,o)a4 + ^-2,2«3 + W-2,o a i + ^0,2^2 = cr ((2|ai | 2 + 2|a 2 | 2 + 2|a 3 | 2 + |a 4 | 2 )a 4 + 2a 3 a!a 2 ) . 

This system (with the derivative terms in y\ = ex\ and y 2 = ex 2 ) was derived in [Tj by using asymptotic 
multi-scale expansions. Higher-order resonances for ui n with larger values of n G N 2 may involve more 
than four components in the coupled-mode equations (|4.6p . and the count of d$ versus n is not available 
in general. 

Remark 3 Additional resonances were considered in p], which correspond to an oblique propagation 
of the resonant Fourier modes, e.g. ei pxi and ea ((p+ 2m i) x i+ 2m 2 x 2) w ithp = —{m\ + m\)/rni ^ Z. These 
resonances can be incorporated in the present analysis by using the transformation U(x) = e^ pXl U{x). 

Remark 4 Existence of two-dimensional (N = 2) gap soliton solutions can not be proved with the 
approach of Section 3 for small values of e when u is close to u n . Indeed, if the solution U (fc) is split into 
a finite number of parts compactly supported near the points of resonances (fci, k 2 ) = \(ni, n 2 ) and the 

remainder part {7o(I&)> then the operator |/c| 2 — w 2 with w 2 = ^j- is not invertible in a neighborhood 
of the circle of the radius \k\ = Since only finitely many parts of the circle are excluded from the 
compact support of Uo(k), the Implicit Function Theorem can not be used to prove existence of the 
map from the finitely many resonance parts of the solution U{k) to the remainder part Uo(k). This 
obstacle has a principal nature as it is related to a generic non-existence of gap solitons in the systems 
without spectral gaps. Indeed, the operator L = —V 2 — eW(x) has no gaps for sufficiently small e in 
the space of two dimensions (N = 2) |17j . 

Remark 5 The conclusions of Proposition Q] and Remark U] can be extended to iV > 3. 
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